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A measure of cluster size heterogeneity ( H ), introduced by Lee et al [Phys. Rev. E 84, 020101 
(2011)] in the context of explosive percolation, was recently applied to random percolation and to 
domains of parallel spins in the Ising and Potts models. It is defined as the average number of 
different domain sizes in a given configuration and a new exponent was introduced to explain its 
scaling with the size of the system. In thermal spin models, however, physical clusters take into 
account the temperature-dependent correlation between neighboring spins and encode the critical 
properties of the phase transition. We here extend the measure of H to these clusters and, moreover, 
present new results for the geometric domains for both d = 2 and 3. We show that the heterogeneity 
associated with geometric domains has a previously unnoticed double peak, thus being able to detect 
both the thermal and percolative transition. An alternative interpretation for the scaling of H that 
does not introduce a new exponent is also proposed. 


I. INTRODUCTION 

Cluster or domain size distributions are commonly 
used in statistical mechanics to unveil geometric prop¬ 
erties and characterize both the equilibrium Q} critical 
behavior and the out of equilibrium dynamics of 

several models. If the system is finite, not all possi¬ 
ble domain sizes are present on a single configuration. 
These sample-to-sample fluctuations disappear and the 
distributions become dense in the infinite size limit or 
after ensemble averages are taken. In order to charac¬ 
terize these fluctuations, a quantity associated with how 
heterogeneously sized the equilibrium domains are, H, 
was recently proposed in the context of explosive per¬ 
colation and then applied to ordinary percolation Q 
and to the Ising jg] and Potts models [?]. Differently 
from the cluster size distribution, the heterogeneity H 
only takes into account whether a given size is present 
in each configuration and gives the number of such dis¬ 
tinct sizes. Although the number of equal-sized clusters 
and their actual size do not enter in the measure of H, 
both entropic and thermal effects are relevant. At high 
temperatures, thermal noise breaks large clusters and the 
fragments are small, and so is the diversity and H. On 
the other hand, for low temperature, the presence of a 
very large, perhaps percolating, cluster leaves small room 
to the other clusters, also decreasing the diversity. 

We here revisit the original formulation for the equi¬ 
librium scaling behavior of H for the Ising model and 
present new results for both geometric and physical clus¬ 
ters Hi- While the former considers all nearest neigh¬ 
bor parallel spins belonging to the same domain, the 
latter only takes the fraction of spins effectively corre¬ 
lated. This difference, although important for systems 
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of interacting spins under thermal noise, obviously does 
not exist for percolation models. For random percola¬ 
tion 0, H has a clear peak and an excellent data col¬ 
lapse is obtained. For the Ising model, no peak has been 
observed for the system sizes considered in Ref. [g] . Al¬ 
though in 2d both kinds of spin domains percolate at the 
same temperature, only the physical ones encode the crit¬ 
ical properties (exponents). It is thus interesting to see 
how their heterogeneity, H p , differs from H g measured 
with the geometric domains. As we show below, even 
for small sizes, H p presents a clear peak and an excellent 
collapse is obtained. For the geometric domains, on the 
other hand, we observe two peaks for sufficiently large 
systems, none of which were observed in Ref. Q, where 
only an abrupt change in the derivative was present at 
the critical temperature T c . The first peak is small, very 
close to T c , while the second one is large, broad and ap¬ 
pears far above the critical temperature. Although the 
very existence of such double peak is interesting enough, 
there are several questions that must be answered. What 
is the origin of such peaks? Do they merge, leaving a sin¬ 
gle peak, as the system size diverges? Do they share the 
same exponents or the percolation exponents play a role? 
What happens in 3d where the critical temperature does 
not coincide with the percolation threshold of geomet¬ 
ric domains? In order to try to answer these questions, 
we studied the equilibrium Ising model on square and 
cubic lattices with linear sizes up to 2560 and 250, re¬ 
spectively. Averages up to 200 samples were performed 
for the smaller systems while larger sizes required fewer 
samples. All measures were obtained after equilibrating 
the system through 500 Swendsen-Wang steps E3- 

The paper is organized as follows. We start Sec. ED by 
reviewing the definition of H and its scaling behavior in 
the context of a thermal system. Then, in Secs. Ill Al and 
III Bl we present the results of extensive simulations for 
the 2d and 3d Ising model, in which we measure H for 
both geometric and physical clusters. We finally discuss, 
in Sec. mu our results and conclude. 
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II. SIZE HETEROGENEITY 


Differently from the susceptibility that corresponds, in 
percolation terms, to the mean cluster size (ID ], and thus 
correlates with both the size and frequency of fluctua¬ 
tions, for the measure of H it is only the presence of a 
given cluster size that is relevant. Both at too low and 
too high temperatures, the distribution of cluster sizes 
is concentrated at small domains, the latter because the 
large noise dissolves the clusters and the former because 
a giant cluster dominates the system leaving almost no 
room for the others. In both cases the size distribution 
has a small tail, domains are rather homogeneous and H 
is not large. We thus expect a maximum of H between 
these two limits, close to the critical temperature because 
the fragmentation of the infinite cluster creates several 
small domains and the size distribution per site, n(s), 
being broad (power law), will increase the probability of 
these clusters having different sizes and thus contribute 
to H. Moreover, as the linear size L increases, more and 
larger clusters are allowed when the average distribution 
is broad enough and H must also increase with system 
size. Notice, however, that having a broad distribution 
does not necessarily imply a large value of H: it depends 
on how many of these different sizes are present on a 
given state or, equivalently, the amount of vacancies in 
the distribution for a particular configuration of a finite 
size system (the average cluster size distribution, on the 
other hand, being an ensemble average, is dense). 

The precise dependence of H on the linear size L can 
be understood with a scaling argument 3,[5j. Let us start 
with the height of the peak close to the critical temper¬ 
ature. Small clusters are present with a high probability 
in a finite size system and we define so, which is a func¬ 
tion of the temperature, as the (average) smallest size 
not present in a configuration. Below so the distribu¬ 
tion is dense, that is, the expected number of clusters 
whose sizes are smaller than so, L d n(s < so), is larger 
than one. On the other hand, above so there are some 
sizes without a realization, L d n(s > So) < 1 and va¬ 
cancies appear in the distribution. Thus, by definition, 
L d n{sg) ~ 1. As we will see, the contributions to H from 
clusters in these two regions are of the same order Q. 


Close to the transition, we write (11) n{s) ~ s T f(z ) 


where z = ts a and t = T/T c — 1 is the reduced temper¬ 
ature. Two critical exponents are introduced. The first 
one, cr, is associated with the extension of the power- 
law tail as the criticality is approached while the Fisher 
exponent r is the exponent of such power-law. For the 
2d cases considered here, we have r g = 379/187 for geo¬ 
metric domains [H, [11], t p = 31/15 for physical clusters 
and T perc = 187/91 for random percolation. The scal¬ 
ing function f{z) approaches a constant for |;z| -C 1 and 
has a fast decay for |z| 3> 1. Close to the transition, 
n(s o) ~ Sq T ~ L~ d , and thus So scales as so ~ L d / T [I[. 
The remaining contribution to H comes from the clus¬ 
ters with s > so- These clusters may be ranked by size, 
starting from the largest SHI- At the criticality, the 


largest cluster size lies in the interval ]si,oo) such that 
the expected number of clusters in this interval is unity. 
That is, the integral of n(s ) from si to oo, multiplied 
by the size of the system is 1. Analogously, the second 
largest cluster lies in the interval [s 2 ,oo) such that the 
integral of n(s) from S 2 to oo, multiplied by the size of 
the system, is now 2. Repeating the process, for the r-th 
cluster in the rank, at the critical temperature, we get 

POO 

r = L d ds s~ T f( 0) ~ L d s].~ T . (1) 


In the appendix we present a more detailed discussion 
of the convergence of this integral. We get that s r ~ 
(L d /r) 1 ^ T ~ 1 ' > and the fractal dimension of these large 
clusters is d/(r — 1) a larger than d/r of sg. As r 
increases, remembering that we are ranking the sizes in 
reverse order, s r — > sg. Since so ~ L d / T , for very large rg 
one has L d l T ~ (L d /rg) 1 ^ T ~ 1 \ hence rg ~ L d / T . Thus, 
the two contributions to H : H ~ So + rg, scale in the 
same way and H ~ L d ^ T Q. Since H and So share the 
same scaling behavior, measuring the first missing cluster 
size So is enough to obtain the critical scaling behavior 
of H. 

The point at which the distribution deviates from the 
power law grows as s* ~ |t| _1 / <T , what defines the expo¬ 
nent a. As pointed out in Ref. [6|, the competing scales 
in this problem can be considered so and s* instead of 
the usual L and £. When s* So, the power law is well 
developed and the distribution for a single configuration 
has many holes (in the opposite limit, s* < so, such dis¬ 
tribution is dense and the system behaves as if its size was 
infinite). In this case, H ~ sg ~ L d l T . Away from the 
transition, s* decreases and the critical region extends 
up to the point in which s* ~ So- Thus, L d s~ T ~ 1 and, 
using the above definition, the size of the critical region 
is |t| ~ L~ da / T = w here = r/ad or, equiva¬ 

lently, zz H = vt /(r — 1). Consistently with this, Lee et 
al 3 considered that the scaling behavior of H was given 
by H(t,L) = L d / T ft.(|f|L 1/,1/H ). While close to the tran¬ 
sition the scaling function h[x) is a constant, away from 
T c the behavior is h(x) ~ x~ x ! a . 


A. 2d 

The definition of H and its scaling were proposed and 
verified in Ref. Q for explosive percolation and, later on, 
also in ordinary percolation Dj and some 2d spin mod¬ 
els @0 for which there are several possible definitions 
of domains. For the Ising model, on which we focus here, 
a (geometric) domain is defined as a connected group of 
nearest neighbors aligned spins (we do not distinguish, 
as in Ref. [6], between up and down spins clusters). We 
show, in Fig. [TJ the geometric domains heterogeneity H g 
as a function of the temperature for several system sizes 
near the 2d Ising critical point. For the small sizes consid¬ 
ered in Ref. 0], only a change in the derivative of H g , very 
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FIG. 1: Size heterogeneity H for the d = 2 Ising model in equilibrium as a function of the temperature for several system 
sizes (top) and the corresponding collapses (bottom) using the reduced temperature t = T/T c — 1. The vertical dotted line 
locates the critical temperature T c = 2/ln(l + %/2). (Left) First (small) peak for the geometric domains, very close to T c . 
The best collapse, keeping v = 1, was obtained with r ~ 2.016(4). (Middle) Second, large and broad peak for the geometric 
domains. The height of the peak grows with an effective r ~ 2.048 that is definitively different from the first peak. A good 
collapse is obtained with v ~ 1.25. The horizontal lines correspond to an independent measure of H at T —» oo using a random 
configuration. (Right) The same but for the physical domains whose exponents are r p = 31/15 and = 31/16. In this case, 
there is a single peak close to T c . Notice the different horizontal scales in all graphs and the logarithmic one in (b). 


close to T c , was observed. By considering larger sizes, in¬ 
terestingly, two well separated peaks develop. The first, 
small one (Fig. |T^), is very close to T c and, as shown by 
the collapse in Fig. [TJi, grows as L d l T . A good collapse 
was obtained with r ~ 2.016(4) (keeping v = 1) and, ac¬ 
cordingly, jz h ~ 1.984. These values are close to the ex¬ 
ponents for the geometric domains in the Ising model |6|, 
r g = 379/187 ~ 2.027 and v n = 379/192 ~ 1.974. On 
a broader range of temperatures (notice the logarithmic 
scale in Fig. [ljo), there is a second peak that, together 
with the first one, was not previously reported. In the 
large temperature limit, H —>• H g { oo) and H g ( oo) can 
be independently measured on a random configuration 
with, on average, half of the spins up and the other half 
down. They are shown as horizontal lines whose height 
slowly increase with L, H g ( oo) ~ logL 11 , and are 
approached as H g (T) — H g ( oo) ~ T -1 . This excess het¬ 
erogeneity at large temperatures for the geometric do¬ 
mains is due to the entropic contribution of uncorrelated 
spins. Because of this size dependence, that persists even 
at very high temperatures where £ —> 0, the scaling of 
this second peak turns to be rather complicated. In¬ 
deed, it does not collapse well with exponents obtained 
for the first peak and, instead, a satisfactory collapse 
was obtained with r ~ 2.048 and v ~ 1.25. These val¬ 
ues are close to the random site percolation exponents, 
r perc = 187/91 and v = 4/3 (thus, = 187/72), but this 


may be just a crossover because of the large separation 
between the two peaks for finite systems (see the 3d case 
below). It seems that even larger systems would be nec¬ 
essary to clarify the nature of the critical exponents in 
this problem. Both peaks, besides increasing in height, 
move towards T c in the thermodynamic limit. Indeed, 
both horizontal scalings in Figs. [TJ:1 and [TJ^ considered 
the known value of T c in d = 2. This is a consequence 
of the fact that, in d = 2, the geometric clusters per- 
colative and the thermal transitions occur at the same 
temperature. 

Since H obviously depends on the nature of the clus¬ 
ters, further insight is gained by comparing with the be¬ 
havior of physical clusters. If we remove from the ge¬ 
ometric clusters the fraction of spins that, albeit paral¬ 
lel, are not effectively correlated, we obtain the so-called 
Coniglio-Klein droplets i m whose physical proper¬ 
ties are also related with the random cluster model |§j. 
Fig-111' shows H p as a function of the temperature for sev¬ 
eral sizes and presents a clear peak close to the transition 
that gets larger as the size of the system increases. No¬ 
tice that the temperatures of all the three peaks discussed 
here are well separated for the system sizes considered. 
Differently from H g ( oo), that is an increasing function 
of L : H p —> 1 when T —> oo since all spins become un¬ 
correlated and form clusters of unit size. The scaling 
proposed in Ref. 0] works very well also for H p , Fig. |Tf, 
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with the Ising thermal exponents (v = 1, r p = 31/15 
and i/ H = 31/16). There is, however, another possible 
point of view for the scaling in this problem. Obviously, 
if \t\L l / vu is a scaling variable, so is (|t|T 1 /'' H ) r /( T_1 ). In¬ 
deed, using |t| r /( T ~ 1 )L 1 / 1 ', collapses equivalent to those 
in Fig. H] are obtained (an example is shown in Fig. [2j, 
with no need to resort to a new exponent vh■ In order 
to collapse the data for H , besides understanding how its 
peak scales with L , one must also know how the width of 
the critical region depends on the system size. Standard 
finite size scaling considers the region around T c in which 
the measured correlation length £ will differ from its infi¬ 
nite size limit behavior. In a finite system, the divergence 
of £ as T —> T c , £ ~ \T — T c \~ u , is constrained by the lin¬ 
ear size L, and the width of the critical region is thus 
\T — T c | ~ L~ v . This interval is, in principle, different 
for other divergent quantities at T c . If we consider H in¬ 
stead of £, since \/v = da/(r — 1) > da /r, the divergence 
of H,H ~ L d / T , is weaker and the critical region, wider. 
Indeed, as discussed in section|TT] \T — T c \ ~ L~( T ~ 1 )l Tl/ _ 



FIG. 2: Collapsed H p for d = 2, analogous to Fig. [TJ, but 
using the alternative scaling with v = 1. 

Fig. [3] shows the fraction of the whole system occupied 
by the largest cluster as a function of the temperature, 
distinguishing whether the largest cluster does or does 
not percolate. At very high temperatures both measures 
coincide since the noise is too high to allow a percolat¬ 
ing cluster. As the temperature decreases, spins become 
more correlated, the clusters larger and these measures, 
still coinciding, steadily increase. They become different 
when some samples percolate at temperatures roughly 
below the second peak of H g . The average fraction size 
of the largest component (top curve) increases almost lin¬ 
early from this point, as T decreases, up to T c where there 
is an abrupt jump and it approaches unity below T c . It is, 
on the other hand, the other measure, the fraction of the 
largest non-percolating cluster that seems to encode the 
relevant information for H and presents a similar double 
peak structure. When some samples start to percolate, 
the largest non-percolating cluster is actually the second 
largest in these samples. As a consequence, its measure 
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FIG. 3: Size of the largest (top) and the largest non¬ 
percolating (bottom) clusters, normalized by N — L d , as a 
function of temperature for L = 640. A cluster is considered 
percolating when it touches all sides of the square lattice. 
Below T c , indicated by a vertical dashed line, there is almost 
always a large, percolating cluster and the available space for 
the second largest cluster is small. For large temperatures, it 
is almost certain that there is no percolating cluster and both 
curves coincide. The insets show snapshots for L = 320 with 
the largest non-percolating geometric cluster in gray at two 
different temperatures (that roughly correspond to the local 
maxima of the bottom curve). At the critical temperature, 
left snapshot, there is a percolating cluster (not shown) and 
the second largest is compact. In red we show, among the 
several physical subclusters in which it divides, the largest 
one. The right snapshot shows that, at higher temperature 
(second maximum), although the largest non-percolating ge¬ 
ometric cluster has almost the same average size as at T c , the 
physical subclusters are much smaller. 


starts decreasing, both measures diverge, leaving a peak 
close to the temperature of the second peak of H. The 
insets of Fig. [3] show, in light gray, one example of the 
largest non-percolating cluster at the temperatures indi¬ 
cated, one close to T c and the other at the point of sep¬ 
aration of the two measures discussed above. These are 
geometric clusters formed by parallel spins. Each one is 
a collection of physical clusters, the largest one is shown 
in red (dark gray). Obviously, as the temperature de¬ 
creases and the spins are more correlated, these clusters 
increase. 


B. 3d 

The heterogeneities can also be measured in the 3d 
case, where the percolative transition of the geometric 
domains does not coincide with the thermal transition, 
T^c < T c ~ 4.51. In 3d, the random site percolation 
threshold for a cubic lattice is 0.312. For all temper¬ 
atures above T c (even in the T —>• oo limit), the frac¬ 
tion of each spin is 0.5, on average, and since both spins 
densities are above the threshold, both percolate. For 
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FIG. 4: (Left) Heterogeneity H for both geometric (left peak, filled symbols) and physical clusters (right peak, empty symbols) 
for the 3d Ising model. The vertical lines are located at T ptsrc ~ 0.95T C (pure geometric percolation) and T c ~ 4.511 (Ising critical 
temperature). Inset: the behavior of H g at the thermal critical point, crossing over to a different declivity. (Right) Collapse 
of both H g (top) and H p (bottom). For H g , t = (T — T perc )/T pBrc , and the best collapse was obtained with T pt>rc ~ 0.95T C , 
r ~ 2.08 and v ~ 0.75. For H p , on the other hand, t = (T — T c )/T c , and an excellent collapse is obtained with 3d Ising model 
exponents: v ~ 0.875, r ~ 2.189. 


temperatures below T c , the system is magnetized and 
the majority spin obviously percolates. The percolation 
transition occurs when, coming from low temperature, 
the minority spin (the one opposite to the system mag¬ 
netization) whose domains are small and isolate, first per¬ 
colate. We obtain, from the best collapse of H (see be¬ 
low) T ~ 0.95T C , slightly larger than the estimate in 
Ref. [H, T perc ~ 0.92T C . 

Fig.[U(left) shows both H s and H p as a function of tem¬ 
perature. Again, H p shows a single peak slightly above 
T c ~ 4.51. On the other hand, as in the 2d case, H g is 
sensitive to both transitions. While there is a large peak 
very close to T perc , at T c there is a crossover to a differ¬ 
ent declivity, as shown in the inset of Fig. 01 While in 
the 2d case two peaks with the same infinite size limit 
were present, in 3d there is a single peak and a crossover, 
each one located in a different temperature for L —> oo. 
The high temperature limit of H g seems to increase even 
slower than in the 2d case, and is barely visible in the 
figure. 

The data collapses in 3d are shown in Fig. 0] (right). 
In the bottom part, the results for H p are very well col¬ 
lapsed, as in the 2d case, with the thermal exponents 
of the 3d Ising model and T c ~ 4.51: v ~ 0.875 and 
r ~ 2.189. In the top figure, for H g , the best collapse 
was obtained with T — 0.95T C , r ~ 2.08 and v ~ 0.75. 
One would expect 17|, in this specific case, the exponents 
of the 3d random site percolation problem, r pero ~ 2.189 
and v pprc ~ 0.875. The reason for this jlj] is that T perc 
being smaller than T c , the correlation length is finite at 
that temperature and clusters at distances larger than 
this length are uncorrelated. The small difference be¬ 
tween this set of exponents and those th at g ive the best 
collapse may be due to finite size effects 18]. 


III. CONCLUSIONS 

The cluster size heterogeneity H is the number of dis¬ 
tinct cluster sizes in a given sample configuration jj] . We 
presented here extensive simulations for the Ising model 
for both geometric (H g ) and physical (H p ) clusters in two 
and three dimensions. While for the latter, H has a single 
peak that diverges at the critical temperature of the Ising 
model, in the former case H presents a signature for both 
the percolative and the thermal transition. Indeed, finite 
2d systems present two quite distinct peaks for H g when 
the size is sufficiently large, but they merge as L —>- oo. 
On the other hand, in 3d there is a peak at T perc and a 
sudden change in the declivity at T c > T perc . Moreover, 
since the peaks of H g and H p do not occur at the same 
temperature in 3d, in the region T perc < T < T c , while 
H p increases, H g decreases (a similar effect occurs in 2d, 
but it disappears in the limit L —>• oo). We also discussed 
the scaling of H and how its behavior correlates with the 
large clusters in the system. Interestingly, although the 
physical clusters collapse is clearly related to the Ising 
thermal exponents, the behavior of the largest peak of 
H g is more subtle. In both 2d and 3d the best collapse 
was obtained with exponents that are close, albeit differ¬ 
ent, from those of the respective random site percolation. 
Although in 3d this is expected (since T perc ^ T C1 the 
correlation length is finite and distant clusters become 
uncorrelated |17|). in 2d both transitions coincide and 
the reason for the similarity with the random site expo¬ 
nents is still an open question. We conjecture that it may 
be just a crossover: for large systems the two peaks get 
closer to each other and converge to T c , the correlation 
length diverges and even clusters far apart are correlated. 
In this way, the exponents should converge to those of the 
first peak. 
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Both in two and three dimensions the heterogeneity 
measured for geometric domains, H g , is very sensitive 
to the thermodynamical transition, even if the critical 
properties are actually encoded in the physical domains. 
The small peak and the declivity change present in 2d 
and 3d, respectively are very close to T c even for very 
modest sizes, and its location does not suffer from strong 
finite size effects. If proven general, H g may provide a 
quite precise method to obtain a preliminary estimate for 
the thermal transition. 

There are, nonetheless, several remaining questions. 
For example, how heterogeneous are the slices of a 3d 
system [18|? Differently from the whole volume behav¬ 
ior, the percolation threshold for the slice again coin¬ 
cides with T c . Another question is whether the above 
behavior of H g occurs in other models as well. We are 
presently studying the Potts model, both for continuous 
and discontinuous transitions, in order to better under¬ 
stand which are the conditions for having a double peak. 
Finally, it is important to perform other geometric mea¬ 
sures in order to better understand the underlying mech¬ 
anism responsible for the behavior of H. 

Appendix: On the convergence of Eq. ® 

When listing the domain sizes present on a particular 
configuration, we notice that all sizes up to a specific size 
appear, while above it several holes are present. After 
averaging over several configurations, we call this size Sq. 

The ranking ro of the cluster size so is 

OO 

r o = L<1 (A- 1 ) 

So 

not necessarily at the critical point. Near it and using 
the scaling relation n(s) ~ s~ T f(z) where z = ts a , one 
could try to transform this discrete sum into an inte¬ 
gral ds s~ T f(z). Unfortunately, this integral does not 
converge inside the critical region |t| ~ ~ L~ da / T . 


However, one can write the n-th derivative of the quoted 
sum 

i fl OO 

- L d J2s n ”- T f (n \z) (A.2) 

SO 

T d i C°° i 

= — / d zz n ~ l -^f^\z). 

(J I ~ 

J Z o 

Now, inside the critical region the bottom limit zq van¬ 
ishes, and the resulting integral converges to a finite (not 
f-dependent) constant for some high enough integer n. 
Therefore, besides additional non-singular terms (pro¬ 
portional to t°, t 1 ... f" _1 ), the leading singular form 
for ?'o reads 

r 0 ~ L d t {T ~ l)/a ~ L d / T . (A.3) 

Notice that the particular mathematical form of the scal¬ 
ing function f(z) does not enter in the argument. The 
above scaling relations are valid anyway, only the propor¬ 
tionality prefactors being dependent on f(z). In Ref. 
the authors assume f(z) to be an exponential function 
and explicitly use this form to derive the above scaling 
relations. However, in 3 dimensions for instance, f(z) 
is very well fitted to a (non centered) Gaussian func¬ 
tion [191 ]. In short, the nice and surprising theory intro¬ 
duced in Ref. Q is correct, but the particular assumption 
adopted by the authors, in order to justify it, is unnec¬ 
essary. 
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